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ABSTRACT 

We present a detailed analysis of the planetary nebula M4-18 (G146.7+07.6) and 
its WClO-type Wolf-Rayet central star, based on high quality optical spectroscopy 
(WHT/UES, INT/IDS, WIYN/DensPak) and imaging (HST/WFPC2). From a non- 
LTE model atmosphere analysis of the stellar spectrum, we derive Tcff=31kK, 
log(M/M0 yr~^)=-6.05, foo=160 km s^^ and abundance number ratios of H/He<0.5, 
C/He=0.60 and O/He=0.10. These parameters are remarkably similar to He 2-113 
([WCIO]). Assuming an identical stellar mass to that determined by De Marco et al. 
for He 2-113, we obtain a distance of 6.8kpc to M4-18 (E'e-v^O.SSmag from nebular 
and stellar techniques). This implies that the planetary nebula of M4-18 has a dy- 
namical age of ^3 100 years, in contrast to >270 years for He 2-113. This is supported 
by the much higher electron density of the latter. These observations may only be 
reconciled with evolutionary predictions if [WC]-type stars exhibit a range in stellar 
masses. 

Photo- ionization modelling of M4-18 is carried out using our stellar WR flux dis- 
tribution, together with blackbody and Kurucz energy distributions obtained from 
Zanstra analyses. We conclude that the ionizing energy distribution from the Wolf- 
Rayet model provides the best consistency with the observed nebular properties, al- 
though discrepancies remain. 

Key words: stars: individual: M4-18; stars: Wolf-Rayet; planetary nebulae: general. 



^ 1 INTRODUCTION 



Low excitation Wolf-Rayet (WR) central stars of planetary 
nebula (PN; denoted [WR] following van der Hucht et al. 
1981) are thought to represent the beginning of hydrogen- 
deficient central star evolution, just following the ejection of 
the PN which occurs at the top of the asymptotic giant 
branch (AGB). However, discrepancies between observed 
abundances and those predicted by theory, led to the pro- 
posal that H-deficient central stars of PN (CSPN) are the 
result of the re-birth of a white dwarf, after a late helium- 
shell pulse (Iben et al. 1983). Following this interpretation, 
the search began for characteristics common to planetary 
nebulae with WC nuclei, that could be used to distinguish 
them from those with H-rich central stars, and thus estab- 
lish that the two classes have followed different evolutions. 
However, it soon emerged that nebulae associated with [WC] 
central stars are indistinguishable from those around H-rich 
CSPN (Corny & Stasinska 1995) . 

To seek a solution to this problem and determine the 
evolution of WR type nuclei, a flurry of empirical and mod- 
elling analyses of WC central stars were carried out by 
several groups (e.g. Leuenhagen, Hamann & JefFery 1996, 
hereafter LHJ). The present study focuses on M4-18 (PN 



C146.7-H07.6, IRAS 04215-^6000), a PN with a [WCIO] cen- 
tral star (following the classification scheme of Crowther, 
De Marco & Barlow 1998). One may reasonably ask why is 
yet another study of a late WC-type CSPN necessary? Re- 
cent stellar (LHJ) and nebular (Surendiranath & Rao 1995, 
hereafter SR) models for M4-18 exist in the literature, to- 
gether with empirical determinations of its PN properties 
(e.g. Coodrich & Dahari 1985). The justification is that pre- 
vious studies considered the star and nebula of M4-18 in 
isolation, making it impossible to relate nebular and stel- 
lar results; additionally a comparison between the nebular 
properties of M4-18 with other PN associated with spectro- 
scopically similar stars suggests that [WCL] stars can follow 
diflterent evolutionary paths; this key aspect has never been 
given due attention^. De Marco, Barlow & Storey (1997) and 
De Marco & Crowther (1998, hereafter DC) carried out a 
rigorous analysis of the central stars and PNe of the [WCIO] 



* Pottasch (1996) remarked that the PNe around M4-18 and 
He 2—113 (another [WCIO] star) are very different, hinting that 
the two central stars cannot have followed the same evolution- 
ary path, but nobody took this argument past Table 10 of his 
contribution. 



© 1999 RAS 



2 0. De Marco & P.A. Crowther 



central stars CPD-56°8032 and He 2-113. We follow these 
methods for M4-18, so that the relative results can be com- 
pared with confidence. 

Our analysis of M4-18 also extends the photo-ionization 
modelling of CPD-56°8032 and He 2-113 by DC using flux 
distributions appropriate for Wolf-Rayet stars. DC identified 
a major discrepancy between the observed and predicted 
nebular properties for these PNe; either the lack of heavy 
element line blanketing in the wind models was to blame, or 
the geometry and high nebular densities of those PNe meant 
that they represented poor probes of the Lyman continuum 
flux of their central stars. As we shall demonstrate, the PN of 
M4-18, with a significantly lower electron density, provides 
a more rigorous test of the theoretical flux distribution. 

In Section ^ we will describe the observations, while in 
Section ^ we will discuss basic observational quantities, in- 
cluding reddening. Section ^ discusses the distance towards 
M4-18, while a quantitative analysis is carried out in Sec- 
tion p. Archive Hubble Space Telescope (HST) images are 
presented in Section ^ while we carry out a nebular abun- 
dance analysis in Section Q In Section ^ photo-ionization 
modelling of the PN is carried out. Finally, we draw our 
conclusions in Section 3. 



2 OBSERVATIONS AND DATA REDUCTION 

We have obtained high quality optical spectroscopy of M4- 
18 using the 4.2 m Wilham Hershell Telescope (WHT), 
the 2.5 m Isaac Newton Telescope (INT) and the 3.5 m 
Wisconsin-Indiana- Yale-NO AO telescope (WIYN). Archival 
International Ultraviolet Explorer (lUE) ultraviolet (UV) 
spectroscopy and HST optical imaging of M4-18 were ob- 
tained from the Uniform Low-Dispersion Archive (Talavera 
1988) at the Rutherford Appleton Laboratory, and the HST 
data archive at the Space Telescope Science Institute, re- 
spectively. The HST Wide Field and Planetary Camera 2 
(WFPC2) data set is discussed further in Section B 



of 40. The data were reduced in a standard manner using 
the IRAF package, with wavelength calibration achieved rela- 
tive to a comparison Th-Ar arc. For flux calibration, the Oke 
(1990) standard G191B2b (DAO) was used. Subsequent data 
reduction was carried out with the dipso package (Howarth 
& Murray 1991). 

Additional high resolution (R~19 000) spectroscopy 
was obtained on November 13, 1998 by D. Sawyer with the 
WIYN 3.5 m telescope and the DensPak fiber array, in the 
ranee 5770-6010A. Three exposures each lasting 600 s were 
obtained with a signal-to-noise ratio of 20. The data were 
reduced with standard iraf routines. Wavelength calibra- 
tion was with respect to a Th-Ar lamp. These spectra al- 
lowed us to measure the radial velocity shifts of the inter- 
stellar Na I D lines. 

Finally, our high resolution INT and WHT observations 
were scaled to the continuum level of the wide-slit INT spec- 
tra, which were obtained during photometric conditions. By 
convolving our observed spectrophotometry with suitable 
synthetic filters (courtesy of J.R. Deacon) we obtained mea- 
surements of both wide-band Johnson photometry (V=14.11 
and B=14.24mag) and narrow-band Smith (1968) photom- 
etry (v=14.14 and fe=14.16mag). Our measurements are in 
reasonable agreement with Shaw & Kaler (1985) who ob- 
tained a (nebula corrected) visual brightness of ~ 14.0 mag. 



2.2 Ultraviolet spectroscopy 

We utilised four (SWP and LWR/LWP) low resolution lUE 
large aperture spectroscopic observations of M4-18 obtained 
between August 1980-August 1991 (LWR8401 was not used 
because it is affected by a cosmic ray hit at the wavelength of 
the C ll] line A2326, Goodrich & Dahari 1985). Although of 
poor quality, these provided constraints on the interstellar 
reddening towards M4-18, and allowed us to obtain esti- 
mates of the nebular carbon abundance. 



2.1 Optical spectroscopy 

M4-18 was observed at the INT, using the Intermediate Dis- 
persion Spectrograph (IDS), together with the 235-mm cam- 
era and a 1024x1024 pixel Tektronix CCD, between July 
17 and 23, 1996. Four settings with the 1200B/Y gratings 
provided complete wavelength coverage in the range 3800- 
68OOA, at a spectral resolution of I.5A using a 1.5 arcsec slit. 
Additional observations with the 300V grating and an 8 arc- 
sec slit (including the entire nebula) provided an absolute 
flux calibration at a spectral resolution of 9.5A. The data 
were reduced in a standard manner using the iraf package, 
with wavelength calibrations achieved using Cu-Ne and Cu- 
Ar arc lamps and absolute flux calibrated data obtained by 
comparing our wide slit observations of M4-18 with obser- 
vations of the Oke (1990) standards B2 iv star BD+33°2642 
and the Op star BD-|-28°4211. 

High spectral resolution [R ~30 000) observations of 
M4-18 were obtained at the WHT using the Utrecht Echelle 
Spectrograph (UES) in service mode on November 27, 1996. 
A 1024x1024 pixel Tektronix CCD provided a complete 
spectral range of 4200-5850A: two exposures each of du- 
ration 2400 sec achieved a continuum signal-to-noise ratio 



3 BASIC OBSERVATIONAL QUANTITIES 

3.1 Radial and nebular expansion velocity 

The radial velocity was measured from Gaussian fits to 
the nebular Balmer lines observed in our optical datasets. 
WHT-UES observations of H/3-7 indicated a mean heliocen- 
tric radial velocity of -51.6±1.0 km s~^ (corresponding to an 
LSR radial velocity of -50.8±1.0km s~^). These results are 
supported by our lower resolution INT observations, which 
indicate a heliocentric radial velocity of -46.1±5.0km s~^. 

We used the FWHM of nebular profiles (H/3, H7 and 
A5755 [N II]) in our WHT-UES observations (FWHM ~6 
kms"'^) to determine the expansion velocity of the PN of 
M4-18, which was revealed to be Wcxp = 19±0.5 km s~^. [The 
value reported on the ESO PN catalogue (Acker et al. 1992) 
of 12 km s~^ derives from an average of [O III] (7.5 km s-i) 
and [N 11] (17.0 km s~^); we note that it is inappropriate to 
average values obtained from lines belonging to a mixture of 
ionization stages, since high ionization lines may originate 
from inner, slower parts of the PN.] 
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Figure 1. (Upper panel): observed ultraviolet and optical spec- 
trophotometry of M4— 18 (units are ergcm"^ from lUE and 
INT, togetlier witli tlieoretical energy distributions reddened by 
Eb-v=0-45 (dashed), 0.55 (dotted), and 0.65 mag (dot-dashed) 
using a standard Galactic extinction curve (Seaton 1979; Howarth 
1983). (Lower panel): As above, except including near and mid- 
IR photometry from SR and IRAS (filled-in squares) , plus a two- 
component blackbody fit (dotted lines) to the IR excess caused 
by dust in the PN. 



Table 1. Balmer fluxes (ergs cm ^ s ^) measured for the nebular 
hydrogen lines of M4-18. 



Line 


WHT-UES 


INT-IDS 


INT-IDS 




narrow slit 


narrow slit 


wide slit 


Ha 




2.69x10"^^ 






6.61x10"" 


5.75x10"" 


8.95x10"" 


H7 


2.30x10"" 







3.2 Interstellar reddening 

Reddening determinations towards M4-18 from the litera- 
ture range from £b-v=0.48 (SR) to £b-v=0.90 (Goodrich 
& Dahari 1985). We derive interstellar reddenings using 
observed Balmer line fluxes, which are listed in Table hi 
Our wide-slit INT H/3 nebular flux is in reasonable agree- 
ment with the determinations of 1.04xl0~^^ ergcm~^s~^ 
by Goodrich & Dahari (1985), and 1.17x10"^^ ergcm"^ s"^ 
by Carrasco, Serrano & Costero (1983, 1984). 

We adopt Ne = 10* cm"^ and Te=10'' K (these are re- 
vised slightly in Section ^), and the Galactic extinction law 
of Howarth (1983) plus hydrogen recombination coefficients 
from Storey & Hummer (1995). Unfortunately Ha is severely 
blended with [Nil] in our wide slit INT observations, so 
we have to rely on narrow-slit observations which imply 
-Bb-v=0.45 mag from Ha-H/J (WHT H/3-H7 observations 
imply i5B-v=0.61). 

We can also derive a reddening towards M4-18 following 
Milne & AUer (1975), using our observed H/3 flux, together 



with the observed 5GHz flux of 22 mjy from Aaquist & 
Kwok (1990). This method implies a higher reddening of 
i5B-v=0.65 mag. 

Therefore, extinctions in the range i5B-v=0.45- 
0.65 mag are implied from the nebula, so we can now test 
which extinction reproduces the observed UV and optical 
stellar flux distribution, making use of our theoretical en- 
ergy distribution from Section ^. In Fig. 0, we compare ob- 
served UV and optical spectrophotometry of M4-18 with our 
theoretical model for interstellar extinctions of iJB-v=0.45, 
0.55 and 0.65 mag. We flnd that 0.55±0.05 mag, provides 
the best match to observations. (The only previous redden- 
ing determination towards M4-18 from spectrophotometry 
was by LHJ who obtained i5B-v=0.70 from a continuum fit 
to exclusively UV data.) We therefore adopt £'b-v=0.55 as 
appropriate to both the nebular and stellar observations of 
M4-18. 

Fig. also includes blackbody fits to the infra-red (IR) 
excess of M4-18, representing warm dust from its PN, based 
on near-IR observations from SR, plus colour corrected 12- 
100/im photometry from IRAS. We find that a two com- 
ponent blackbody fit is necessary, indicating the presence of 
hot (~1800K) and warm (~300K) dust. DC found somewhat 
cooler dust properties for CPD-56°8032 and He 2-113. 



4 THE DISTANCE TO M4-18 

Distance estimates to M4-18 range from 1-7 kpc (SR; Cahn, 
Kaler & Stanghellini 1991). De Marco et al. (1997) obtained 
estimates of the distances to CPD-56°8032 and He 2-113 
based on the LSR radial velocities of the interstellar Na I D 
lines, or following the assumption that the entire bolomet- 
ric luminosity is re-radiated in the IR by cool circumstellar 
dust. The electron density of the PN of M4-18 is more than a 
factor of ten lower than that of CPD-56°8032 and He 2-113 
(Section with IRAS fluxes about 30 times lower, sug- 
gesting that the stellar luminosity is unlikely to be entirely 
re-radiated in the infrared. 

In Fig. ^ we present the WIYN/DensPak observations 
in the spectral region 5880-5900A, showing the Na i D lines 
superimposed on a stellar emission blend (C 11 M5, C ill M20 
and He 11 A5897, De Marco et al., 1997). To determine the 
radial velocities of the individual Na I D line components, 
we used a model that calculates absorption proflles for a 
variety of interstellar clouds with Gaussian line-of-sight ve- 
locity distributions (implemented into the ISCALC routine in 
the DIPSO package). To rectify the Na I D lines we 'snipped' 
the D2 absorption and neglected the spectrum longward of 
~5892A, thereby approximating the broad emission feature 
with a single Gaussian (dashed line in the upper panel of 
Fig. |). 

The three-cloud component fit to the Na I D2 line is 
presented in Fig. ^ (lower panel), with the corresponding pa- 
rameters listed in Table H Although we were unable to carry 
out the same procedure for Dl, it is clear from Fig. |^ that 
the positions of the Dl absorption line components match 
those of the D2 components. 

If the -60 km s~^ radial velocity component is of inter- 
stellar origin, the galactic radial velocity map of Brand & 
Blitz (1993) indicates that M4-18 lies beyond 10 kpc. How- 
ever, at this distance, M4-18 would lie 1.2 kpc above the 



© 1999 RAS, MNRAS 000, §-|l| 



4 



0. De Marco & P.A. Crowther 




5884 5888 5892 
Wavelength (A) 



5896 




-160 



-80 
Velocity (km/s) 



Figure 2. Upper panel: The Na I D lines in the WIYN spectrum 
of M4— 18. The dashed line shows a one-Gaussian fit to the main 
stellar emission feature. Lower panel: The rectified D2 spectrum 
in the LSR rest frame, fitted by a three-cloud model (dashed line). 
Dl is also shown. The vertical line marks the LSR radial velocity 
of the PN. 



Table 2. Parameters of the three-cloud model shown in Fig.^ 
The velocity dispersion parameter b is defined in Howarth & 
Phillips (1986), while A'^ is Na I the column density. The radial 
velocity (RV) components are in the LSR rest frame, while D is 
the distance implied by the Brand & Blitz (1993) galactic rotation 
curve for material near the galactic plane. 



Cloud no. 


h 


Log(Af/cm ^) 


RV 


D 




(km s^-"-) 




(km s-i) 


(kpc) 


1 


13.0 


13.0 


-60.0±5.0 




2 


20.0 


12.8 


-37.0±5.0 


>1.5 


3 


12.0 


13.6 


-10.0±5.0 


>1.0 



galactic plane, where neutral sodium is likely to be scarce 
(the scale height of neutral hydrogen was determined to be 
144±80 pc by ShuU & Van Steenberg 1985), Moreover, the 
Brand & Blitz (1993) radial velocity curve was obtained 
from measurements of H ll regions with a scale height of 
only 67 pc. In the direction of M4-18, it can be considered 
reasonably accurate up to a distance from the galactic plane 
of a few times this value («200 pc; several measurements of 
H II regions at ~400 pc above the galactic plane exist for this 
line-of-sight). Interstellar material in the direction of M4-18, 
therefore extends to a distance of ^1.5 kpc, corresponding 
to a radial velocity of -30 km s"'^. Na l D line components 
with radial velocity more negative than this value, cannot 
be interstellar. Considering the uncertainty on the Galactic 
rotation curve, we identify the observed -37 km s"'^ compo- 
nent to be of interstellar origin. Consequently, Na I D lines 
simply argue for the distance towards M4-18 being greater 
than 1.5 kpc. 



We attribute the -60 km s~ component to nebular 
material. Although we may expect M4-18 not to have a 
neutral envelope since it is optically thin (see Section 
Na I was observed in other PNe by Dinerstein, Sneden & 
Uglum (1995), some of which are also optically thin. If Na I 
is present in the PN shell, its expansion velocity can be ob- 
tained from the difference between the LSR radial velocity 
of the PN (-50.8 km s~^, Section ^ and the radial velocity 
shift of the Na I D line component at -60 km s~^. This re- 
sults in 9 km s^^, significantly lower than the ionized shell 
expansion velocity (19 km s~^). Neutral shells moving con- 
siderably slower than the ionized envelopes were detected 
by Dinerstein et al. (1995) in two other PNe (SwSt 1 and 
IC3568). 

Therefore, we follow the alternative approach of DC 
for CPD-56°8032 and He 2-113 based on assuming a core 
mass of 0.62 Mq (obtained by X.-W. Liu [priv. comm.] 
for five LMC WR CSPN), and applying the helium burn- 
ing post-AGB evolutionary tracks of Vassiliadis & Wood 

(1994) to obtain a stellar luminosity of 5250 L©. Utilising 
the observed «-band magnitude for M4-18 from our spec- 
trophotometry, together with the interstellar extinction of 
£'B-v=0.55=b0.05 m ag a nd bolometric correction from our 
WR model (Section 5.2) we obtain a distance of «6.8kpc. 
This is in excellent agreement with recent statistical dis- 
tances for M4-18, which range from 6.7-7.1 kpc, as sum- 
marised by Zhang (1995) and with the distance implied by 
the -37 km s~^ Na I D line component. 

Of course, the masses of all WC central stars need not 
be identical. From Vassiliadis & Wood (1994) and Blocker 

(1995) , the possible range in luminosity for a CSPN is 
~2 500-16 000 L0, corresponding to a possible range in dis- 
tance of 4.8-12.0 kpc. Consequently, all distances lower than 
this minimum are excluded. For example, SR adopted a dis- 
tance of 1 kpc, based on a relationship between the dust 
temperature and the PN radius (Pottasch et al. 1984). Us- 
ing our bolometric correction, an unrealistically low stellar 
luminosity of HOLq would be implied for this distance. (Al- 
ternatively, an unrealistically high bolometric correction of 
—6 mag would be required assuming the luminosity from SR 
namely 2 5001/© ) . 

In summary, we adopt a distance of 6.8 kpc to M4-18, 
assuming a mass and luminosity identical to that obtained 
by DC for He 2-113. However, ahhough M4-18 and He 2- 
113 are spectroscopically very similar (Section |^, this as- 
sumption could be in error. The existence of massive WC 
stars and CSPN with very similar spectra, such as He 2-99 
([WC9]) and HD 164270 (WC9) (see Mendez et al. 1991), is 
a clear indication that spectral similarities do not necessar- 
ily indicate similarities in mass. This estimate, on the other 
hand is in agreement with the limits imposed by stellar evo- 
lution. 



5 STELLAR ANALYSIS 

In this Section we will present our stellar spectrum of M4- 
18, addressing the important question of whether stellar hy- 
drogen is present. We will also obtain stellar parameters 
using a sophisticated stellar atmosphere model appropriate 
for Wolf-Rayet stars. 
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Figure 3. Comparison between rectified spectra of M4-18 (WHT-UES, WIYN-DensPak for A>5845A) and He 2-113 (DC), demon- 
strating the striking similarity of their appearances. 



5.1 The stellar spectrum — is hydrogen present? 

The stellar spectrum of M4-18 is dominated by emission 
lines of helium (He l-ll), carbon (C ll-iv) and oxygen (O ll- 
lll). Our high resolution spectroscopy of M4-18 is presented 
in Fig. ^ along with a spectrum of the [WCIO] central star 
He 2-113 (from De Marco et al. 1997). It is apparent that 
their spectral morphologies are remarkably similar, includ- 
ing line widths, suggesting comparable outflow velocities. 
The nebular lines of the M4-18 PN are much stronger than 
in He 2-113. 

The upper panel of Fig. ^ shows the high resolution 
rectified UES spectrum, in the rest frame of the H/3 line. This 
feature appears to broaden at its base; additionally, a blue- 
shifted absorption (whose minimum intensity is registered 
at 185 km s~^ to the blue of the line's rest wavelength), is 
apparent from the rectified spectrum in Fig. ^ (upper panel), 
suggesting that stellar hydrogen is present in the wind. LHJ 
determined a H/He ratio of <10, 0.64 and <0.5 by number, 
for the [WCIO] central stars M4-18, He 2-113 and CPD- 
56°8032, respectively. De Marco et al. (1997) noticed the 
same behaviour in the Balmer lines of CPD-56°8032 and 
He 2-113. However, by comparing these profiles with nebular 
[O l] lines at 6300 and 6363A, they demonstrated that the 
origin of the broad pedestal under the Ha and H/3 profiles 
is more likely to be due to an irregular nebular geometry. 

For M4-18, nebular [N ii] A5755 and H7 are also present 
in the high resolution UES spectrum, although these are too 
weak for a comparison to be made. However, we are able to 
comment on the origin of the P Cygni profile, which cannot 



be produced by peculiar nebular material. The lower panel 
of Fig. ^ compares the P Cygni absorption feature in the rest 
frame of He 11 4859. ISA (transition 8-4), together with the 
adjacent Pickering members A4541.46 (9-4) and A5411.37 
(7-4). We find that Hen is responsible for this feature. 

To summarise, we suspect that H/3 emission wings are 
solely due to high velocity components in the nebular spec- 
trum, as is the case for the other two [WCIO] central stars 
He 2-113 and CPD-56°8032. (De Marco et al. 1997; Sahai, 
Wotten & Clegg 1993). To verify this, one would need to 
observe other nebular lines at high spectral resolution. 



5.2 Spectroscopic analysis 

Hillier (1987, 1990) theoretical model atmospheres are used 
in our spectroscopic analysis. These calculations employ an 
iterative technique to solve the transfer equation in the co- 
moving frame subject to statistical and radiative equilibrium 
in an expanding, spherically-symmetric, homogeneous and 
steady-state atmosphere. The stellar radius (R*) is defined 
as the inner boundary of the model atmosphere and is lo- 
cated at a Rosseland optical depth of 20. The temperature 
stratification is determined from the assumption of radiative 
equilibrium with the temperature parameter (T,) defined by 
the usual Stefan-Boltzmann relation. Similarly, the effective 
temperature CTeft) relates to the radius (-R2/3) a-t which the 
Rosseland optical depth equals 2/3. The spectral synthesis 
proceeds by fitting line profiles to observed diagnostic Hel- 
11, C ll-iv, O ll-lll lines, together with the absolute visual 
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Table 3. Derived stellar parameters of M4-18, including a comparison with Leuenhagen et al. (1996, LHJ) scaled to our assumed 
luminosity, and demonstrating its very close similarity with He 2-113 (from De Marco & Crowther 1998, DC). We include the predicted 
Hi and Hel continuum ionizing fluxes (Qo, Ql)- 



Star 


Study 


V 


d 








log L, 


log M 




H/He 


C/He 


O/He 


log Qo 


log Ql 


A/v 






mag 


kpc 


mag 


kK 


Rq 


Lq 


Me yr~' 


kms^^ 








s-i 


s-i 


mag 


M4-18 


This work 


14.1 


6.8 


0.55 


31 


2.4 


3.72 


-6.0 


160 


<0.5 


0.60 


0.10 


47.2 


36.7 


-1.8 




LHJ 


13.3 


3.5 


0.70 


31 


2.3 


3.67 


-6.0 


350 


<10 


0.40 


0.05 






-1.6 


He 2-113 


DC 


11.9 


1.2 


1.00 


31 
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3.72 


-6.1 


160 


0.0 
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36.7 


-1.8 
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Figure 4. (upper panel) Base of the rectified H/3 profile in the 
WHT-UES dataset of M4-18, revealing the presence of (i) broad 
line wings and (ii) a P Cygni component; (lower panel) Compar- 
ison of rectified Hen Pickering series profiles, indicating that the 
H/3 P Cygni component is due to Hen A4859 (8-4). 



magnitude. In the absence of high resolution ultraviolet ob- 
servations we measured a stellar wind terminal velocity from 
optical Hel P Cygni profiles, resulting in 160±15kms~^. 

Fig. P presents a comparison of selected observed 
(WHT-UES) line profiles (solid) with our synthetic spec- 
tra (dotted) for M4-18. Overall, the quality of these and 
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Figure 5. A comparison between selected theoretical profiles for 
M4-18 (dotted lines) with WHT-UES observations (solid lines). 
Other profile fits are comparable to those shown for He 2-113 in 
DC. 



other profile fits, covering a wide range in excitation and ion- 
ization, is comparable with that achieved for other [WCIO] 
stars by DC, where a more comprehensive explanation of 
the successes and failures of the model is presented. Over- 
all the fits successfully reproduce the strength and shape of 
many optical emission line profiles. In particular the helium 
spectrum is well fitted, although it is clear from the poor 
fit to the width of He ll A4686 (Fig. || - top right panel) 
that the adopted l3=l velocity law is not ideal for this star. 
P Cygni absorption components are poorly predicted (e.g. 
C IV AA5801-12), further supporting an inadequate velocity 
structure. 

Table ^ presents a summary of our derived stellar pa- 
rameters. Overall, as refiected in their spectra, the stellar 
properties and chemistries derived for M4-18 and He 2- 
113 are almost identical, except for a marginally stronger 
stellar wind and higher carbon content in the former. For 
reasons discussed by DC we consider that derived param- 
eters are not critically dependent on our neglect of heavy 
element line blanketing. Although we favour a non-stellar 
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origin for the pedestal at the base of the H/3 hne, we have 
determined a strict upper hmit for the stellar hydrogen 
abundance from fitting the wings of the broad H/3 feature, 
yielding H/He<0.5, by number. Our upper limit, imposes 
a tighter constraint on the hydrogen abundance than the 
analysis of LHJ; this is solely due to the higher spectral res- 
olution of our observations (10 kms~^ versus ~150 kms~^). 

LHJ included M4-18 in their quantitative study of 
[WCL] CSPNe. Although their study used an independent 
code, the same assumptions of spherical symmetry and ho- 
mogeneity were made. Table |^ also compares our derived 
stellar parameters with those from LHJ. Overall agreement 
is very good, although LHJ adopt a higher terminal velocity, 
again due to the lower resolution of their data set. LHJ de- 
rive somewhat lower carbon and oxygen contents for M4-18. 

Note that the distance to M4-18 adopted by LHJ is in- 
consistent with our value, despite deriving an almost iden- 
tical stellar luminosity. This is because LHJ obtained an 
extinction of i?B-v=0.7mag on the basis of a spectral fit to 
exclusively UV data, inconsistent with our stellar and neb- 
ular data sets. Their distance estimate of 3.5 kpc (following 
Cudworth 1974) then led to a visual continuum magnitude 
of v=13.3 mag, which our spectrophotometry does not sup- 
port (Section 



6 THE M4-18 NEBULA 

We now discuss the PN associated with M4-18, including 
archive HST/WFPC2 narrow-Ha (F656N) images, obtained 
on April 18, 1996 by R. Sahai. From previous ground based 
observations, the PN of M4-18 is known to be extremely 
small. Radii from optical, mid-IR and radio data sets sug- 
gest ^2" (Shaw 1985; Zijlstra, Pottasch & Bignell 1989). 
Although the HST/WFPC2 Ha image has recently been 
presented by Sahai & Trauger (1998), we show it in Fig. ^ 
in connection with further measurements of its size which 
are relevant to our analysis. We will abstain however from 
discussing its morphology, which is described in detail by 
Sahai & Trauger (1998). Fig. ^ shows the nebula both as a 
grey-scale and a contour map, indicating a clear elongated 
ring which peaks in intensity at a semi-major axis of 1.28" 
(PA=0) and a semi-minor axis of 0.8" (PA=90). The ring 
is embedded in a larger, less eccentric shell. From an az- 
imuthal average of the PN, we determine a mean radius of 
1.85" ±0.10". 

At a distance of 6.8 kpc, the radius of M4-18 corre- 
sponds to a physical radius of 0.06 pc (10"* AU). Using the 
observed nebular expansion velocity we obtain a dynamical 
age of rdyn~3 000 years. Note that this dynamical age should 
be reasonable since as we will later show, the PN is optically 
thin (Section |^). 



7 NEBULAR ABUNDANCE ANALYSIS 

Spectroscopy of the M4-18 PN has previously been carried 
out by Sabbadin (1980), Goodrich & Dahari (1985) and SR. 
The latter groups combined low resolution optical and UV 
spectroscopy to determine PN physical conditions and abun- 
dances, while SR also applied their own nebular modelling 



code. We shall now re-derive the electron density and tem- 
perature for the PN of M4-18, together with elemental abun- 
dances. 

T.l Nebular line fluxes 

Our analysis of the PN associated with M4-18 is based on 
our INT spectrum, in which stellar and nebular features 
are blended. (Since the nebular diameter is ~3.7", ground- 
based long-slit spectroscopy do not allow off-star nebular 
extraction) . We present observed and de-reddened line fluxes 
in Table ^ 

The [O ll] doublet at 3726 and 3729A may be contami- 
nated by a component of the stellar O ll M2 triplet A3727.3. 
This was judged to be negligible since the un-blended com- 
ponents of this triplet (at 3712.7 and 3749. 5A) were weak. 
However, this doublet was located close to the end of our 
spectral range, where the sensitivity of the instrument drops 
sharply. We therefore assign a 30% error to our measure- 
ments. Inspection of our UES data set reveals that both 
the [Olll] lines at 5007 and 4959A are present^, though ex- 
tremely weak. 

Following the approach discussed in De Marco et al. 
(1997) for the C ill] A1909 line, we consider this to be entirely 
of stellar origin. Instead, we rely on C ll] A2326 for the de- 
termination of the nebular carbon abundance. Its observed 
equivalent width is 15±4A from our dereddened lUE spec- 
trum, of which ~4A is of stellar origin according to our stel- 
lar analysis (Section Since the continuum flux at 2326A 
is in good agreement with our reddened theoretical model 
distribution (Fig. hi) , we are able to obtain a nebular A2326 
flux of 5.34xl0~^*erg cm~^ s~\ 

For sulphur, only [S ll] AA6713,31 were available. The 
[S III] A3722 line is extremely weak (and blended with [O ll] 
AA3726,29), while the AA9069,9532 lines were outside our ob- 
served spectral range. SR derived the abundance of S^^/H"*", 
although they derived fluxes for AA9069,9532. from indirect 
techniques. Finally, for nitrogen we used lines at 5755 and 
6548A, excluding [N ll] A6584 which is blended with the 
stellar C ll line at 6578A. 

7.2 Nebular temperatures, densities and 
abundances 

We can now proceed to obtaining estimates of the nebular 
conditions and chemistry. From the usual diagnostic dia- 
gram relating Hq/[Sii] to [Sii] A6717/A6731 (see Sabbadin, 
Minello & Bianchini 1977) we flnd that M4-18 falls in the 
photo-ionization dominated region, rather than being ion- 
ized by shock waves. Therefore the usual nebular diagnostic 
techniques may be applied here. 

We have constructed a nebular diagnostic diagram for 
M4-18 using [Nil], [Sll] and [On] line ratios. These are gen- 
erated using the RATIO program, written by I.D. Howarth 

t The [Olll] line at 4959A is heavily blended with C II stellar 
components (i.e. the four components of the dielectronic Multiplet 
25), but as we have fitted this spectral region in the spectra of 
CPD-56°8032 and He 2-113 which did not have any [O III] (De 
Marco, Storey & Barlow 1998) it was straightforward to recognise 
its presence. 
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Table 4. Observed (F) and dereddened (I; _Ba_^/=0.55 mag) 
nebular line intensities for M4-18 from our INT/IDS spec- 
trum (*: WHT/UES). The dereddened H/3 flux is 3.62x10-^2 



erg cm 



Errors are indicated. 



Ion 


A 


F 




100 X I 


Errors 




(A) 


(ergs cm" 


-2 s-i) 


mm 


% 


[O II] 


3726.0 


5.69x10" 


13 


158 


30 


[O II] 


3728.8 


3.01x10" 


13 


83.7 


30 


H/3 


4861.3 


5.75x10" 


13 


100 


10 


[O III] 


5006.8* 


<3.6xl0 


-15 


0.06 


50 


[Nil] 


5754.6 


1.11x10" 


14 


1.35 


10 


[O I] 


6300.3 


2.15x10" 


14 


2.2 


10 


[O I] 


6363.8 


9.73x10" 


15 


1.0 


10 


[Nil] 


6548.0 


4.44x10" 


13 


43.1 


10 


H« 


6562.8 


2.69x10" 


12 


260 


10 


[N II] 


6583.4 


1.54x10" 


12 


147 


20 


[S II] 


6716.5 


7.19x10" 


14 


6.7 


10 


[S II] 


6730.8 


1.41x10" 


13 


13.0 


20 
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Figure 7. Nebular Diagnostic diagram for M4-18. 



and S. Adams which solves the equations of statistical equi- 
librium allowing a determination of as a function of for 
each ratio. From Fig. ^ we determine A''^ =io^-^='='''2 cm"^ 
and Te=8 200±200 K for M4-18. The electron density is 
adopted as a compromise between the [O ii] and [S ii] dou- 
blets. Our results are in reasonable agreement with Sabbadin 
(1980) and SR , though Goodrich & Dahari (1985) derived 
a somewhat lower electron temperature (re=5600±1600 K, 

We now utilise our derived nebular parameters to ob- 
tain estimates of elemental abundances, which are listed in 
Table ^. For oxygen, nitrogen and carbon we believe the 
abundances derived from the singly ionized ions to be rep- 



resentative of the whole ionized populations since [O ill] 
A5007 and C lll] A1909 are extremely weak or absent. In 
the case of oxygen, from the upper limit of A5007, we derive 
O2+/H+<4.4xl0-'', corresponding to O2+/O+<1.0x 10"^ 
Equally for sulphur, S^/H"*" should be representative 
of S/H since SR obtained S^+/S+~0.1. Our C/0 ratio is 
much lower than found by De Marco et al. (1997) for CPD- 
56°8032 and He 2-113, which suggests that the carbon abun- 
dance in the PNe of WC central stars is not always higher 
than average. 

Overall, agreement with the abundances obtained by SR 
is good, except for oxygen (they derived log(0/H) + 12=8.1) 
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Table 5. Nebular abundances for M4-18 (on the usual scale 
where logN(H)=12.0). The C/O number ratio is on a linear scale. 
The [Neil] line flux is taken from Aitken & Roche (1982). 



Ratio 


log(X/H)+12.0 


Diagnostic Lines 


C 


9.08±0.03 


C ll] A2326 


N 


7.60±0.04 


[N II] A5755,A6548 


O 


8.62±0.01 


[O ii] AA3726,29 


Nc 


6.73±0.05 


[Ne ii] 12.8 Aim 


S 


6.16±0.06 


[S II] AA6713,31 


C/O 


2.9±0.6 





since their de-reddened [O ll] AA3727,29 flux was almost 10 
times lower than our measurement. Abundances derived by 
Goodrich & Dahari (1985) are also in reasonable agree- 
ment once the different reddening and electron properties 
are taken into account. 



8 NEBULAR MODELLING 

In this Section we derive Zanstra temperatures for M4-18 
using blackbody and Kurucz (1991) model atmosphere flux 
distributions, which are compared to WR flux distribution 
resulting from the analysis of the CSPN (Section The 
nebular radial density profile is derived from the HST Ha 
image surface brightness distribution following the formal- 
ism of Harrington & Feibelman (1983). We will finally carry 
out photo-ionization modelling of the PN of M4-18 using 
all available atmosphere fiux distributions and we will fi- 
nally compare the modeled parameters with those derived 
empirically in Section |^. 



8.1 Zanstra temperatures 

We derive a fiux of hydrogen-ionizing photons correspond- 
ing to log Qo=46.81 s"^ from our H/3 fiux (Table |l), de- 
reddened by iJB-v=0.55 at a distance of 6.8 kpc. We derive 
a blackbody H I Zanstra temperature of 29 500±500 K, while 
Kurucz ATLAS9 (Kurucz 1991) model atmospheres yield 
an effective temperature of 33 000±500 K, using log g—4:.0 
models (the minimum gravity tabulated at the required ef- 
fective temperature). The inferred bolometric luminosities 
are 4000 Lq and 5600 Lq, respectively. For comparison, 
Goodrich & Dahari (1985) derived a temperature of 22 000 K 
using the Stoy (1933) energy-balance method. SR derived a 
blackbody H I Zanstra temperature of 26 OOOK. (They also 
fitted a blackbody to the optical and UV energy distribution 
implying 23 000 K.) 

Our WR stellar model for M4-18 predicts log Qo=47.23 
s~^. This is 2.6 times greater than that inferred from the 
nebular H/3 flux. In contrast, DC showed how the WR model 
atmosphere ionizing flux of He 2-113 was 85 times larger 
than that predicted by H/3, suggesting that for high density 
PN, gas-dust competition might play a role. Our results for 
M4-18 confirm their suspicion. 



1 - 




12 3 

R(arcsec) 

Figure 8. H/3 surface brightness profile (in units of 
erg cm~^ s"'^ arcsec"^) derived from the HST Ha image of M4- 
18 from an azimuthal average centred on the star (solid line; the 
central star's continuum is responsible for the intensity peaking 
near the centre of the nebula). The dashed line is the intrinsic 
surface brightness distribution predicted by our photo-ionization 
model. 



8.2 Nebular surface brightness and density 
distributions 

We have obtained the radial density distribution for M4-18, 
necessary for photo-ionization modelling, from the archival 
HST/WFPC2 image (Section §). For this PN the contribu- 
tion from nebular continuum emission is negligible. 

The azimuthally averaged surface brightness distribu- 
tion in the Ha image was derived using an algorithm imple- 
mentation in the SURFPHOT package of MIDAS. The centre 
of the image was chosen to coincide with the position of 
the central star. The integrated flux within a radius r was 
then determined. The central star was contained within a 
radius of 5 pixels and the contribution of its continuum flux 
was subtracted from the nebular flux. The nebular flux, in- 
tegrated over the entire PN, was then normalised to the 
de-reddened H/3 flux, differentiated with respect to r and 
divided by 27rr to obtain the required azimuthally averaged 
H/3 surface brightness distribution. The result is plotted in 
Fig. ^ (solid line). 

For optically thin PNe, the formalism of Harrington & 
Feibelman (1983) allows one to derive a density distribution 
from a surface brightness distribution (see Liu et al. 1995 for 
an implementation of this method) . Photo-ionization models 
based on density distributions derived with this formalism, 
however, tend to give an average density which is too low 
to reproduce density-sensitive line ratios such as the [O ll] 
A3726/A3729 ratio. This deficiency is overcome by introduc- 
ing a filling factor, e, such that within any small region of 
material a fraction e is flUed with gas, while the complemen- 
tary fraction (1-e) is empty. In our case all density sensitive 
line ratios were well reproduced with a flUing factor of unity 
(i.e. a completely filled nebula). 

The density distribution which reproduces the H/3 sur- 
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Figure 9. (Upper panel): Electron density distribution for M4- 
18 obtained from our analysis of the HST/WFPC2 image; (Lower 
panel): Electron temperature resulting from our photo-ionization 
modelling with the WR flux distribution (assuming a distance of 
6.8kpc). 



face brightness distribution, absolute H/3 value and density 
sensitive line ratios is shown in Fig. M (upper panel). To 
reproduce the surface brightness shown in Fig. N we used 
A'' = 12 and if = 1 in equation A2 of Harrington & Feibel- 
man (1983). The upper panel of Fig.^ shows that M4-18 has 
a central cavity surrounded by a thick shell whose density 
declines fairly symmetrically as the radius increases. 



8.3 The photo-ionization modelling 

The photo-ionization code used to model M4-18 (Harring- 
ton et al. 1982) assumes that the nebula can be represented 
by a hollow, spherical shell which is ionized and heated solely 
by the radiation of the central star. This shell is sampled by 
a network of grid points at each of which the density is 
supplied. We used 60 grid points and the elect ron density 
derived from the H/3 surface brightness (Section 8.2). 

Once a distance is adopted, and the de-reddened V mag- 
nitude of the central star established, the effective temper- 
ature of the central star (from the WR stellar analyis or 
derived by a Zanstra analysis as described in Section i.l) 
determines all the other stellar parameters. These param- 



eters are listed in the first three rows of Table ^. The in- 
ner nebular radius was fixed at the value derived from the 
HST/WFPC2 observations, scaled to the adopted distance. 
Elemental abundances are initially set at empirical values, 
and later adjusted to fit the observed line fluxes if necessary. 
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Figure 10. Comparison of Lyman ionizing flux distribution ob- 
tained from our stellar analysis of M4— 18 (Wolf-Rayet, solid), with 
black-body (dashed) and Kurucz (1991, dotted) distributions ob- 
tained from our nebular Zanstra analysis. Clearly the WR distri- 
bution differs dramatically from the other distributions shortward 
of the He I edge at A504, especially at the On edge at A353. 



Final photo-ionization models for each flux distribution 
are compared to observation in Table ^. Overall, very good 
agreement is reached between observations and model calcu- 
lations for all diagnostic line ratios (O ll, N ll and S ll) as well 
as absolute line fluxes. This is an indication of the success of 
our density distribution, recalling the discrepancy between 
different diagnostics in Fig. |^ 

It is apparent that the WR and Kurucz (1991) stel- 
lar flux distributions perform better than the blackbody. 
Overall, the WR stellar atmosphere is most appropriate to 
reproduce the line intensities, in spite of the other distribu- 
tions resulting from Zanstra nebular analyses. In Fig. ^, we 
present the predicted electron temperature as a function of 
radius for our WR model. Hydrogen is completely ionized 
throughout, as is appropriate for an optically thin PN, with 
a predicted nebular mass of O.OSMq. 

A major failure of each model flux distribution is the 
overestimate of the sulphur ionization balance, and conse- 
quent underestimate of the [S ll] 6717, 3lA line strengths. 
(The ratio S^+/S+ was estimated to be 0.1 by SR, although 
a value of <1 is probably more secure.) The poor repro- 
duction of the sulfur ionization balance can be appreciated 
by inspecting the shape of the ionizing flux distributions in 
Fig. |lC|. Metal line blanketing in our WR model atmospheres 
may help to solve the sulphur problem by reducing the hard 
UV flux (for the effect of blanketing in late-type WR stars 
see Crowther, Bohannan & Pasquali 1998). 

In the case of the oxygen ionization balance, only 
the WR model predicts the appropriate ratio, with Ku- 
rucz (1991) and blackbody atmospheres overestimating the 
strength of [O ill] A5007. This is fllustrated in Fig. ^ where 
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Table 6. Comparison of results from photo-ionization models 
for M4-18 with observations (I(H/3)=100), using WR non-LTE, 
blackbody and Kurucz (1991) plain— parallel energy distributions. 
Elemental abundances are by number. Filling factors of e=l are 
assumed throughout. 



Parameter Empirical 


WR 


Blackbody 


Kurucz 






31.0 


29.5 


33.0 


log(L*/L0) 




3.72 


3.60 


3.75 


log Qo(s"^) 


46.81 


47.23 


46.81 


46.81 


A^e(r)(cm~^) 


See Figure 


y (upper panel) 
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4.0 
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5.0 


4.1 


4.1 


4.0 


o/ n X iu 


1.4 


1.4 


1.4 


1.4 


[O ii]A3726/A3729 


1.88 


1.82 


1.85 


1.80 


[O iij AU 1 -l-U/ AU 1 Ol 


0.52 


0.60 


0.60 


0.62 


[N ii]A5755/A6548 


0.032 


0.027 


0.029 


0.032 


Cm] A1909 





0.0 


11 


2.6 


Cii] A2326 


148 


128 


96 


197 


[On] A3726 


158 


155 


154 


213 


[On] A3729 


84 


85 


83 


118 


[Oiii] A5007 





0.0 


39 


1.0 


[Nil] A5755 


1.4 


1.13 


1.14 


1.35 


[O i] A6300 


2.2 


0.049 


0.067 


0.23 


[O i] A6363 


1.0 


0.016 


0.022 


0.075 


[Nil] A6548 


43 


42 


41 


42 


[Sii] A6717 


6.7 


0.64 


0.87 


1.29 


[Sii] A6731 


13 


1.73 


1.85 


2.08 


[Neil] 12.8;im 


32 


32 


32 


32 


S2+/S+ 


<1 


12 


15 


17 


02+/0+ <1.0xlO"^ 


~0 


77x10^3 


1.6x10" 


Te(N+){K) 


8200 


7900 


8000 


8400 


r(H I) 




1.1 


1.6 


4.2 



only the WR energy distribution has negligible zero flux be- 
low the He I edge at A504. 

Our study supersedes the photo-ionization modelling of 
M4-18 by SR. They adopted a distance of 0.9kpc, and a 
22 000 K blackbody, implying a luminosity of 134L0 for the 
central star. As discussed in Section ^, these quantities are 
inconsistent with theoretical predictions for CSPNe. Note 
also that SR also adopted a constant electron density of 
6620 cm~^, implying Te=6 600 K (their empirical value was 
8 500 K). 



9 DISCUSSION 

We have derived consistent stellar and nebular parameters 
for the [WCIO] central star M4-18 and its PN. Although 
stellar and nebular analyses have previously been carried 
out (LHJ; SR) such studies revealed inconsistent parameters 
(e.g. LHJ disregarded its visual magnitude while SR adopted 
an unrealistically low luminosity). 

While the stellar parameters of M4-18 are almost iden- 
tical to those of He 2-113, their PNe differ greatly. For M4- 
18, a dynamical age of 3 100 years is implied for a distance 
of 6.8 Kpc, while the well determined distance to He 2-113 



Table 7. Summary of stellar and nebular parameters for M4- 
18, compared with those for its spectroscopic twin He 2-113 (De 
Marco et al. 1997; DC). Evolutionary time-scales (tcvoi) are from 
Blocker (1995) and Vassiliadis & Wood (1994). 



D Ne Hcxp Radius r^yn i* M t^voI 

kpc cm""^ km s"-*^ lO^'^pc yr Lq Mq yr 









M4- 


-18 










4.8 


3300 


19 




43 


2200 


2600 


0.52 


8 000 


6.8 


3300 


19 




61 


3100 


5250 


0.62 


500 


12 


3300 


19 


He 2- 


111 

-113 


5600 


16000 


0.84 


150 


1.2 


63 000 


21 




6 


>270 


5250 


0.62 


500 



(De Marco et al. 1997), implies a dynamical age for this PN 
of ~270 yr. Since He 2-113 is optically thick, its dynami- 
cal age is probably only a lower limit, although there can 
be little doubt that the PN of He 2-113 is younger, when 
we consider the respective optical depths and electron den- 
sities. How can such spectroscopically similar CSPN have 
such different PN ages? This difference is hard to reconcile, 
if M4-18 and He 2-113 have similar masses. 

Were we to relax our assumption that all [WC]-type 
stars have an identical mass, distances in the range 4.2- 
12 kpc are consistent with our spectrophotometry and 
the luminosity from theoretical He-burning tracks. Table ^ 
shows the properties of M4-18 for this range of distances 
compared to those of He 2-113 (De Marco et al. 1997; DC). 
We include the ages predicted by the evolutionary tracks of 
Blocker (1995) and Vassihadis & Wood (1994). (Note that 
these agree well where comparisons may be carried out). 
Ages are predicted to be extremely sensitive to core mass. 
From Figures 15 and 16 of Blocker (1995), a helium-burning 
central star of mass 0.52 Mq reaches 30 000 K in 8 000 years, 
while a star of 0.84 Mq takes only 150 years. 

Comparisons between dynamical and evolutionary ages 
should be treated qualitatively (both dynamical and evolu- 
tionary ages of the PNe are subject to large uncertainties). 
Nevertheless, assuming a lower mass for M4-18 than He 2- 
113 resolves the differences in their PNe. In addition, it re- 
solves the discrepancy between the dynamical and evolution- 
ary ages of M4-18. For a mass of '^O.55M0, the dynamical 
and evolutionary time-scales are both ~2500 yr (correspond- 
ing to a distance of 5.5 kpcQ). We note that for He 2-113 
the comparison between dynamical and evolutionary time- 
scales is reasonable if we take into consideration that the 
dynamical time-scale is merely a lower limit. 

If the masses of M4-18 and He 2-113 are different, then 
masses of WR central stars may span the whole range of 
possible post AGB masses, opening the possibility that WR 
central stars follow a variety of evolutionary paths. It has 
been suggested (e.g. Rao 1987; SR) that [WCIO] central 
stars are approaching the AGB in a born-again type evolu- 
tion (Iben et al. 1983). This claim was based on the appar- 
ent inverse correlation between PNe diameters and CSPN 



t At a distance of 5.5 kpc, the luminosity, radius and mass- 
loss rate for M4-18 would be revised to log (L»/L0)~3.45, 
R,=1.8Rq and log(M /Mq yr~l)=-6.25, with other stellar prop- 
erties unchanged. 
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temperatures. Incorporating the results from DC, this is no 
longer established. Consequently, there is no evidence that 
[WCIO] CSPN are approaching the AGB. However, our com- 
parison provides evidence that WC central stars are not a 
homogeneous group, with a range of masses and individual 
evolutionary patterns. It remains to be confirmed whether 
WC-type stars can result from a 'born-again' evolution, al- 
though the potentially related stars Abell 30 and Abell 78 
are associated with well established born-again PNe ( Jacoby 
& Ford 1983). 
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